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ABSTRACT 

In the context of theories of gravity modified to account for the observed dynamics of galactic 
systems without the need to invoke the existence of dark matter, a prediction often appears 
regarding low acceleration systems: wherever a falls below ao one should expect a transition from 
the classical to the modified gravity regime. This modified gravity regime will be characterised 
by equilibrium velocities which become independent of distance, and which scale with the fourth 
root of the total baryonic mass, V 4 cx M. The two above conditions are the well known flat 
rotation curves and Tully-Fisher relations of the galactic regime. Recently however, a similar 
phenomenology has been hinted at, at the outskirts of Galactic globular clusters, precisely in 
the region where a < oq. Radial profiles of the projected velocity dispersion have been observed 
to stop decreasing along Keplerian expectations, and to level off at constant values beyond the 
radii where a < ao. We have constructed gravitational equilibrium dynamical models for a 
number of globular clusters for which the above gravitational anomaly has been reported, using 
a modified Newtonian force law which yields equilibrium velocities equivalent to MOND. We find 
models can be easily constructed having an inner Newtonian region and an outer modified gravity 
regime, which reproduce all observational constraints, surface brightness profiles, total masses 
and line of sight velocity dispersion profiles. Through the use of detailed single stellar population 
models tuned individually to each of the globular clusters in question, we derive estimates of the 
total masses for these systems. Interestingly, we find that the asymptotic values of the velocity 
dispersion profiles are consistent with scaling with the fourth root of the total masses, as expected 
under modified gravity scenarios. 

Subject headings: gravitation — stellar dynamics — stars: kinematics — globular clusters: general 



1. Introduction 

In the context of modified gravity theories, 
where the dynamical measurements usually as- 
cribed to the presence of dark matter are inter- 
preted as evidence for a change in the form of grav- 
ity, stellar populations in the Galactic halo have 
proven to be critical test grounds. The relative 
proximity of many such systems makes them ac- 
cessible to detailed photometric and spectroscopic 
kincmatical observations, while the total absence 
of gas in the cases of many dwarf spheroidal satel- 
lites (dSph) and globular clusters, implies these 
objects present a relatively clean test, as only the 
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measured stellar component must be responsible 
for the measured dynamics, in the absence of any 
dark matter. 

Galactic dwarf spheroidal satellites appear un- 
der classical gravity as the most heavily dark 
matter dominated systems known, and have been 
modelled extensively under MOND e.g. Sanchez- 
Salcedo et al. (2006), Angus (2008), Hernandez 
et al. (2010), McGaugh & Wolf (2010), Kroupa 
et al. (2010) for some recent examples, or equiv- 
alently, thinking of MOND as a modified Newto- 
nian force law, rather than a change to Newton's 
second law, by Mendoza et al. (2011). Sanchez- 
Salcedo & Hernandez (2007) also studied the prob- 
lem of the tidal limiting of Galactic globular clus- 
ters and local dSph galaxies comparatively under 
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MOND and dark matter, finding both hypothesis 
to be consistent with the data, given the present 
level of observational errors. 

One generic prediction of modified gravity the- 
ories designed not to require the hypothesis of 
dark matter, is that as accelerations fall below 
a « 1.2 x 10~ 10 m/s 2 , a transition should occur 
away from Newtonian gravity (e.g. MOND in Mil- 
grom 1983, TeVcS in Bekenstein 2004, QUMOND, 
BIMOND in Zhao & Famaey 2010, modified New- 
tonian force law in Mendoza et al. 2011). This 
modihed regime will be characterised by equilib- 
rium velocities for test particles orbiting within 
spherical mass distributions which become con- 
stant with distance, and which scale with mass 
as V A oc M. 

In this sense, the recent studies by Scarpa et 
al. (2007) and Scarpa et al. (2011) appear par- 
ticularly exciting, as they find precisely the above 
mentioned transition in a number of Galactic glob- 
ular clusters, precisely beyond the radius where 
accelerations fall below a . 

The problem of Galactic globular clusters under 
modified gravity theories has been treated before, 
e.g Baumgardt et al. (2005) study the mean val- 
ues of the velocity dispersion expected under New- 
tonian gravity and MOND comparatively, finding 
larger values are expected in MOND, for globular 
clusters in the outer halo of the Milky Way. Mof- 
fat & Toth (2008) however, obtain expected values 
for the mean velocity dispersion of GCs compat- 
ible with Newtonian expectations, and projected 
velocity dispersion profiles which fall slower than 
Newtonian expectations, being compatible with 
the measurements of Scarpa et al. (2007). Haghi 
et al. (2009) use a MOND N-body code to account 
for the effects of the external field effect of MOND 
in the problem, and find that the profiles of the 
projected velocity dispersion of Galactic globular 
clusters can serve as a test in constraining mod- 
ified gravity theories. Using an analytical treat- 
ment of the problem under MOND, ignoring the 
external field effect of MOND, Sollima & Nipoti 

(2010) construct the MOND equivalent of analytic 
King profiles, where one limits the extent of the 
modelled distributions through the Galactic tides, 
to get models where the velocity dispersion pro- 
files fall to zero. Also very recently, Haghi et al. 

(2011) use again a numerical MOND N-body code 
to show that under MOND mean velocity disper- 



sions of Galactic GCs, are expected to be higher 
than Newtonian. 

In this paper we construct fully self-consistent 
dynamical models for Galactic globular clusters 
using the modified Newtonian force law formula- 
tion of Mendoza et al. (2011), to explore the con- 
sistency of that approach. The modelled density 
profiles are calibrated to match observed surface 
brightness profiles, projected velocity dispersion 
radial profiles, and total masses, inferred through 
a careful single stellar population modelling of the 
globular clusters in question, taking into account 
metalicity and age of the relevant stellar popula- 
tions. We pay particular attention to reproduc- 
ing the flat outer regions of the velocity dispersion 
profiles, as recently measured by Drukier et al. 
(1998), Scarpa & Falomo (2010) and Scarpa et al. 
(2011). 

We find that fully self-consistent equilibrium 
models can be constructed to match all observa- 
tional constraints available on Galactic globular 
clusters, including the outward flattening of the 
projected velocity dispersion profiles. To within 
uncertainties, the same happens under Newtonian 
gravity (e.g. Drukier et al. 2007, Lane et al. 2010, 
Kiipper et al. 2010), where the kinematics in ques- 
tion are interpreted as evidence of tidal heating by 
the overall Galactic gravitational field, occurring 
even in the absence of prominent tidal tails, or 
contamination of unbound stars. It is interesting 
that we further find that the masses and asymp- 
totic velocity dispersions of the studied globular 
clusters are consistent with a scaling <r oc M. 
This last point strengthens the interpretation of 
the observed dynamics in the outskirts of Galactic 
GCs as evidence for modified gravity in general. 

In section (2) we derive the model through 
which equilibrium profiles for spherically symmet- 
ric stellar populations are derived, under the modi- 
fied Newtonian force law of Mendoza et al. (2011). 
In section (3) we show that such models can be 
easily obtained to satisfy all observed parameters 
for a sample of 8 recently observed Galactic glob- 
ular clusters, all showing clearly a flattening of 
their projected velocity dispersion profiles at large 
radius. Our conclusions appear in section (5). 
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2. Non isothermal gravitational equilib- 
rium models 

We shall model globular clusters as populations 
of self gravitating stars in spherically symmetric 
equilibrium configurations, (e.g. Sollima & Nipoti 
2010) under a modified Newtonian gravitational 
force law for test particles at a distance r from 
the centres of spherically symmetric mass distri- 
butions M(r): 



f(x) = a x 



1 



1 - X™- 1 



(1) 



In the above x — lu/r, where Im — {GM{r) / clq) 1 / 2 . 
We see that when a >> a , x » 1 and 
f(x) — > a x 2 , one recovers Newton's gravity 
force law, while when a << ciq, x << 1 and 
f(x) — > a x 7 where an equivalent MOND force 
law / = (GM^ao) 1 / 2 /r is obtained. The in- 
dex n mediates the abruptness of the transition 
between the two regimes. Notice that for a soft 
n = 2 transition function one recovers the classical 
MOND p function of Bekenstein (2004). 

This was shown in Mendoza et al. (2011) 
to yield generalised isothermal gravitational equi- 
librium configurations with characteristic radii, 
masses and velocity dispersions, r g , M, a, which 
smoothly evolve from the classical virial equilib- 
rium of M — <r 2 r g /G to the observed tilt in the 
fundamental plane of elliptical galaxies, to the 
a- 4 = (MGao) scaling of the galactic Tully-Fisher 
relation, as one goes from i » 1 to i ~ 1 to 
x « 1. Also, consistency with solar system ob- 
servations was found there to constrain the tran- 
sition to be fairly abrupt, requiring n > 4. In 
Mendoza et al. (2011) it is also proven that a suf- 
ficient condition for Newton's theorems for spheri- 
cally symmetric mass distributions to hold, for any 
modified force law, is only that / can be written 
as a function exclusively of the variable x. It is 
this last which we will be using in what follows, to 
construct equilibrium models for globular clusters. 
Recently, Bernal et al. (2011) showed this mod- 
ified force law model to be the low velocity limit 
of a formal generalisation to GR of the f(R) type, 
providing a theoretical basis for the model used. 

The equation of hydrostatic equilibrium for a 
polytropic equation of state P = Kp 1 is 



2 d P 



dr 



-v<p. 



(2) 



Since p — (Airr 2 )^ 1 dM (r) / 'dr , going to locally 
Maxwellian conditions 7 = 1 and K = <r 2 (r), the 
preceding equation can be written as: 



dM(r)\ _1 d 2 M{r) 2 



dr 



dr 2 



= —clqx 



(3) 

where a(r) is the isotropic Maxwellian velocity dis- 
persion for the population of stars, which is al- 
lowed to vary with radius, as observed in globular 
clusters, e.g. Sollima & Nipoti (2010). The above 
is a generalisation of the treatment presented in 
Hernandez ct al. (2010), which we used in the 
modelling of dSph galaxies, which are charac- 
terised by flat velocity dispersion profiles, obtain- 
ing mass models consistent with observed velocity 
dispersions, half mass radii and total masses, in 
the absence of dark matter. Locally Maxwellian 
models of this type with radially varying volumet- 
ric velocity dispersions can be found in e.g. Ibata 
et al. (2011), where conditions are further gener- 
alised to the inclusion of a varying radial orbital 
anisotropy, this last, for simplicity, we take as zero. 

As an illustrative example we can take the limit 
f(x) — aox 2 , and recover —GM(r) jr 2 for the 
right hand side of equation (3), the Newtonian ex- 
pression appearing for a » do, or equivalently 
x « Im- If one then imposes isothermal condi- 
tions cr(r) = a and looks for a power law solution 
for M(r) = M (r/r ) m , we get: 



a 



m 



GM 



(4) 



and hence m = 1, the standard isothermal halo, 
M(r) = 2<r 2 r/G, having a constant centrifugal 
equilibrium velocity v 2 = 2a 2 and infinite extent. 
At the other limit, a << a , x » Im, equation 
(3) yields: 



[GM a ]V2 ^ r 
Jo 



m/2 



(5) 



In this limit m = 0, we obtain M(r) — M and 
v 2 = So 2 = (GMoao) 1 / 2 , the expected Tully- 
Fisher scaling of the circular equilibrium velocity 
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with the fourth root of the mass, with rotation 
velocities which remain flat even after the mass 
distribution has converged, rigorously isothermal 
halos are naturally limited in extent. It is interest- 
ing to note that in this limit the scaling between 
the circular rotation velocity and the velocity dis- 
persion is only slightly modified as compared to 
the Newtonian case, the proportionality constant 
changes from 2 to 3, for the squares of the veloc- 
ities. In astrophysical units, this low acceleration 
limit for an isothermal halo in gravitational equi- 
librium yields: 



cr = 0.2 



0.35 



Mo 

Mr 



Mo 



1/4 



km/ s, 



1/4 



km/s, 



(6) 



(7) 



for the velocity dispersion and the centrifugal equi- 
librium velocities of a halo of total mass Mo, re- 
spectively. 

Also, it must be borne in mind that Galactic 
globular clusters are not in isolation, but orbit 
within the Galactic environment. This leads to a 
radius beyond which the tidal forces of the Milky 
Way unbound the outer stars of a globular cluster, 
which depends on the mass and orbital radius of 
the GC in question. The above tidal radius can be 
calculated to first order for a point mass orbiting 
another point mass by equating the internal grav- 
itational force of the GC to the derivative of the 
Galactic gravitational force, we briefly recall the 
first order tidal limit calculation of: 



dF ext {R) 



dR 



Ar = Ft, 



(8) 



which leads to the tidal density stability condi- 
tion of p s > ~p for the density of a satellite of ex- 
tent Ar and mass M s orbiting at a distance Rq 
from the centre of a spherical mass distribution 
M(R) having an average matter density p inter- 
nal to i?o resulting in a gravitational force F ext (i?), 
where the internal gravitational force under New- 
ton is F int — — GM 2 / (Ar) 2 , under the assump- 
tion Ar << R. The equivalent calculation un- 
der the force law given by the a << ao limit of 
F = -(GMa ) 1/2 /R is given by: 



(GM(R)a )^ 2 (GM s a ) 1/2 
Ar = 



R 2 



leading to: 



Ar 
R 



Ps> "FT IP 



(9) 



(10) 



as the equivalent of the classical tidal density cri- 
terion. In terms only of masses and radii we obtain 
the tidal radius as: 



»( Ms V /4 

n = R {mJr)) 



(11) 



Since the spatial extent of globular clusters will 
always be much smaller that their Galactocentric 
radii, equation(lO) shows that under the modified 
gravitational force law explored here, GCs will be 
much more robust to tides than under Newtonian 
gravity. As an example, for a 3 x 10 5 M Q GC orbit- 
ing lOkpc from a 1.6 x 10 11 M Q mass, correspond- 
ing under eq.(7) to the 220km/ s of the Galactic ro- 
tation curve, the tidal radius comes to 370pc. This 
is relevant as it shows that treating GCs as isolated 
systems is a self-consistent assumption, within the 
approach of the modified force law being consid- 
ered. Recent examples of dynamical modelling of 
Galactic GCs treated as isolated systems can be 
found in e.g. Sollima & Nipoti (2010) and Haghi 
et al. (2011). 

Taking initial conditions M(r) — > 0, dM(r)/dr — 
4irr 2 po, when r — > 0, a constant central density 
Po, we can solve the full second-order differen- 
tial equation (|3]) for M(r) through a numerical 
finite differences method, once a model for cr(r) is 
adopted. 

Solving equation ([3]) yields the volumetric pro- 
files for the density and mass, p(r), M(r). p(r) 
is then projected along one dimension to obtain 
a projected surface density mass profile, E(i?). 
In all that follows we shall use r for radial dis- 
tances in 3D, and R for a projected radial coordi- 
nate on the plane of the sky. S(i?) can be com- 
pared to observed surface brightness profiles once 
a mass-to- light ratio is assumed. Given the ap- 
pearance of mass segregation processes in the dy- 
namically evolved stellar populations of Galactic 
globular clusters, volumetrically, the mass-to-light 
ratio of a real GC will never be strictly constant. 
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However, given the ages of the stellar populations 
involved, only a narrow range of stellar masses 
remain, and further, observed surface brightness 
profiles are projections on the plane of the sky of 
the volumetric density profiles. Thus, observed 
surface brightness measurements towards the cen- 
tral regions of GCs imply an integration across 
the entire foreground and background of any ob- 
served cluster. This in turn means that the regions 
where mass segregation is strongest, the central 
ones, contribute to the projected surface bright- 
ness profiles only after being averaged out over 
lines of sight traversing entirely the cluster in ques- 
tion. Even after this averaging process, the high 
densities found towards the centre still surely im- 
ply M/L values which strictly must be a function 
of radius, even for the surface brightness profiles 
I(R). Still, for the reasons described above, it is 
frequent to compare modelled mass surface 

density profiles directly to observed projected sur- 
face brightness profiles through the use of constant 
M/L values, e.g. Jordi et al. (2009), Haghi et al. 
(2011). In what follows we will also perform such 
comparisons assuming constant M/L values. 

Finally, a projected velocity dispersion profile 
<j p (R) is constructed through a volume density 
weighted projection of the volumetric velocity dis- 
persion profile of a model, <r(r). Therefore, a pro- 
posed model for er(r) does not directly give a pro- 
jected a p (R) profile, which the model only yields 
after having solved for the detailed density struc- 
ture, and the subsequent mass weighted projection 
of the proposed cr(r). 

Observations provide only the density weighted 
projected profiles of a(r), integrated along the line 
of sight, <jp(R), not the volumetric a(r) profiles 
which we require in equation ([3]). We shall there- 
fore adopt a parametric form for the volumetric 
velocity dispersion profile, and adjust the param- 
eters to obtain a match to the observed globular 
cluster properties. For this we use: 

a(r) = o-icxp + a o- (12) 

In the above ao is given directly by the obser- 
vations as the asymptotic value of the measured 
projected velocity dispersion profile for each glob- 
ular cluster, as at large R radii, projection effects 
tend to zero and a p (R) — > a(r). This leaves us 
with three model parameters to determine: po , g\ 
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Fig. 1. — The figure shows the volumetric density 
profile for a sample model. The vertical line gives 
the point where x = 1. Notice that the asymptotic 
p(r) profile at large radii is steeper than the r~ 2 
of the Newtonian case, resulting in a finite total 
radius and total mass for the configuration. 

and r CT , which we fit to match the observed pro- 
jected velocity dispersion profiles, as well as cen- 
tral projected velocity dispersion values for each 
observed globular cluster treated, and comparing 
the resulting model S(i?) profiles to the observed 
surface brightness profiles of any given GC, un- 
der the requirement that the M/L ratios used be 
consistent with detailed inferences from the stel- 
lar evolutionary and IMF studies of McLaughlin 
& van der Marel (2005). 

We end this section presenting in figures (1) 
and (2), a sample model. In figure (1) we show 
the volumetric density profile, which is qualita- 
tively similar to a cored isothermal profile, with 
the difference that the asymptotic p(r) profile at 
large radii is steeper than the r~ 2 of the Newto- 
nian case. This results in finite total masses and 
finite half-mass radii even for the asymptotically 
flat cr(r) volumetric profiles we adopt, in contrast 
with the situation in classical gravity, where in- 
finitely extended mass profiles would appear. This 
finite profiles are also what appears under the type 
of modified gravity laws we are treating, even for 
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Fig. 2. — The figure shows the volumetric mass 
profile for the same sample model of fig.(l). The 
vertical line gives the point where x = 1. The 
convergence to a total mass is evident, despite not 
having introduced any externally determined tidal 
radius. 

rigorously isothermal cr(r) = cr equilibrium con- 
figurations, as already pointed out in Hernandez 
et al. (2010) and Mendoza et al. (2011) and as 
shown in the developments following eq. (5). 

Figure (2) gives the corresponding volumetric 
radial mass profile for the same sample model, in- 
dicating again with a vertical line the threshold 
where x = 1. We see that only 20% of the total 
model mass lies beyond x — 1. This fraction is 
distributed over a much larger area having much 
smaller projected mass surface densities than the 
central regions interior to the x = 1 threshold, the 
ones which are much more easily observed, and 
over which Newtonian gravity accurately holds. 

These two figures illustrate the physical con- 
sistency of the model, a positive isothermal volu- 
metric velocity dispersion is assumed, integration 
of eq.(3) then yields a volumetric density profile 
essentially consistent with what would appear un- 
der Newtonian dynamics in the region interior to 
x = 1, where the force law converges precisely 
to the standard expression. At large radii how- 
ever, the density profile increasingly steepens and 



naturally reaches p{r) = at a well defined total 
radius, as is evident from the convergence seen in 
figure (2). Thus, the distribution function is neces- 
sarily positive throughout the modelled structure, 
and goes to zero at a well defined outer radius. In 
the following section we give best fit models for 8 
recently observed Galactic globular clusters. 

3. Modelling observed globular clusters 

We now present best fit models constructed as 
described above, optimised to match observed pro- 
jected surface brightness profiles, under M/L val- 
ues consistent with the ranges given by McLaugh- 
lin & van der Marel (2005) for each individual GC 
modelled, and projected velocity dispersion mea- 
surements for a series of Galactic globular clusters. 
The data for the projected velocity dispersion pro- 
files we take from Drukier et al. (1998), Scarpa 
et al. (2004), Scarpa et al. (2007a), Scarpa et 
al. (2007b), Scarpa & Falomo (2010) and Scarpa 
et al. (2011), and the surface brightness profiles 
from Trager et al. (1995). Total masses we in- 
fer from the integration of the observed surface 
brightness profiles, assuming again M/L values 
consistent with the ranges given by McLaughlin 
& van der Marel (2005) for each individual GC 
modelled, detailed single stellar population models 
tuned to the inferred ages and metallicities of each 
of the clusters we model, and checked for consis- 
tency through the synthetic HR diagram construc- 
tion procedure described in Hernandez & Valls- 
Gabaud (2008). 

For this last, the relevant ages and metalici- 
ties we take from the literature: Salaris & Weiss 
(2002); Kraft & Ivans (2003) for NGC 288 and 
NGC 6341, Salaris & Weiss (2002); Rutledge et al. 
(1997) for NGC 6171, Salaris & Weiss (2002) and 
McNamara et al. (2004); Kraft & Ivans (2003) for 
NGC 7078 and Samus et al. (1995); Harris (1996) 
for NGC 7099, respectively. In obtaining inferred 
total masses, confidence intervals were assigned by 
considering the full range of M/L values given by 
McLaughlin & van der Marel (2005) for each in- 
dividual GC, for the full range of plausible stellar 
models considered in that study, for each of the 
GCs we treat. The final results for the mass to 
light ratios appear in the table at the end of this 
section. 

The final resulting density weighted projected 
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Fig. 3. — The figure compares the resulting model projected surface brightness profiles and resulting model 
projected velocity dispersion profiles to the corresponding observed quantities for the first four globular 
clusters studied. In all cases the assumed M/L values where within the ranges given by McLaughlin & van 
der Marel (2005) for each individual GC through detailed stellar population modelling. The vertical line 
gives the point where x = 1 and the modified force law used shifts from the classical Newtonian form to the 
MONDian character of eq.(5). 
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Fig. 4. — The figure compares the resulting model projected surface brightness profiles and resulting model 
projected velocity dispersion profiles to the corresponding observed quantities for the last four globular 
clusters studied. In all cases the assumed M/L values where within the ranges given by McLaughlin & van 
der Marel (2005) for each individual GC through detailed stellar population modelling. The vertical line 
gives the point where x = 1 and the modified force law used shifts from the classical Newtonian form to the 
MONDian character of eq.(5). 
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velocity dispersion profiles, a p (R), together with 
the resulting surface brightness profiles for the 8 
GCs treated appear in figures (3) and (4), where 
the vertical line indicates the point where a = a^. 
The models accurately fit the observed projected 
surface brightness profiles of the clusters, the cen- 
tral value of the measured projected velocity dis- 
persion, as well as the observed profiles for a p (R). 

We see five cases, NGC 288, NGC 1851, NGC 
5139, NGC 6341 and NGC 7099 where the mod- 
els very accurately reproduce all light distribution 
and velocity dispersion profiles, these last in fact 
agree with the observations in all cases, to within 
reported uncertainties. For NGC 1904 and NGC 
6171 a very slight mismatch appears between the 
models and the observations of the surface bright- 
ness profiles at values of R ~ a few pc, and for 
NGC 7078, over the outer regions, precisely be- 
yond x = 1, the model falls below the observed 
surface brightness profile. These two last points 
could be evidence of the failure of some of the as- 
sumptions used in these cases, probably the pres- 
ence of velocity anisotropy in the dynamics of the 
stars in question. Alternatively, a slight tuning 
of the assumed volumetric M/L ratios with ra- 
dius, within entirely plausible ranges, would im- 
prove the fits. Rather than introduce a further 
degree of freedom, we prefer to show that highly 
adequate fits are easily obtained, for the simplest 
M/L = constant assumption. 

We see also the steepening of the surface density 
brightness profile towards the edge of the clusters, 
particularly for NGC 288 and NGC 7099. Un- 
der Newtonian gravity this feature would be in- 
terpreted as a tidal radius, while under the modi- 
fied force law scheme treated here, it is a natural 
consequence of the change in gravitational regime, 
leading to finite matter distributions when a < do 
and a(r) tends to a constant cf., eq.(5). 

A model with a degree of radial orbital anisotropy 
which varies with radius is an entirely plausible 
alternative, implying the introduction of a further 
free function which allows to more accurately re- 
produce the surface brightness profiles. We have 
also chosen not to complicate the model with the 
inclusion of a /3(r) profile, and preferred to show 
that under the simplest isotropic construction, 
the modified force law we test is capable of ad- 
equately reproducing all the observations in five 
of the cases studied, although slight mismatches 
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Fig. 5. — Dependence of resulting model mass on 
the value of the force law index n. The ratio of the 
model mass to the inferred globular cluster mass 
is shown for two representative clusters, with the 
shaded region giving the range allowed by direct 
stellar population studies of the clusters shown, 
NGC 288 and NGC 6341. This illustrates the need 
for a sharp transition in the force law used. 

appear in three other cases. In particular, for 
NGC 7078, our model requires more mass at large 
radii, but without an increase in a(r) in that re- 
gion. This can be trivially accommodated by an 
increase of kinetic energy at large radius which 
should not be evident in the observed velocity dis- 
persion, through the inclusion of radial motion, 
a degree of orbital radial anisotropy appearing 
beyond 20 pc. 

Obtaining a good match required taking a high 
exponent of n = 10 in equation (3), i.e., a rel- 
atively sharp transition for the generalised force 
law, which however, remains continuous and dif- 
ferentiable at all points, by construction. This is 
in agreement with observed upper limits to devi- 
ations from Newtonian dynamics at the solar sys- 
tem, models for the Galactic rotation curve, and 
equilibrium models for the local dSphs, as shown 
in Mendoza et al. (2011) to apply for the modi- 
fied force law tested here, provided n > 4. Essen- 
tially this implies that the transition between the 
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Newtonian and the "MONDian" regimes must be 
fairly abrupt. In fact, given the form of equation 
(3), the final model becomes independent of n for 
any value of n higher than 10. 

We now illustrate the dependence of the mod- 
els obtained to the assumed force law index n. 
Figure (5) shows the ratio between the resulting 
model mass, M, to the observed cluster mass, M c , 
as a function of the value of n used, for best fit 
models for two illustrative cases, NGC 288 and 
NGC 6341, upper and lower curves at the right 
end, respectively. The shaded region gives the al- 
lowed range of globular cluster mass for a par- 
ticular choice of stellar population parameters in 
McLaughlin & van der Marel (2005), if the sys- 
tematic uncertainties were included, the shaded 
region would extend, but only upwards, by a fac- 
tor of about 2. We see that obtaining a total model 
mass in accordance with direct stellar population 
studies requires taking a value of n > 5 for NGC 
288, and n > 8.5 for NGC 6341. The particular 
minimum value of n required for the other clusters 
in the sample varies somewhat, although most re- 
quire values n > 8 — 9 to reach the minimum in- 
ferred masses. As this parameter can not be ex- 
pected to change from cluster to cluster, we have 
taken n = 10 in all cases. 

Notice that this value results in a equiva- 
lent MOND fx transition function which is fairly 
abrupt. For comparison, most MOND transition 
functions proposed imply a softer transition from 
the Newtonian to the MONDian regimes, e.g. the 
\i function of Bekenstein (2004) correspond ex- 
actly to the force law treated here, with n = 2. 
A similar result has recently been presented by 
Qasem (2011), who analyses laser lunar ranging 
data to constrain the allowed departures from 
Newtonian gravity at the scale studied, and finds 
the data rule out most proposed MOND transi- 
tion functions, while our proposal of Mendoza et 
al. (2011) remains consistent with the test per- 
formed. This completely independent study sup- 
ports the constraints of figure (5), implying that 
the transition from the Newtonian regime to the 
MONDian one is probably steeper than commonly 
thought. 

Figure (6) now gives the radial profiles of a p (R) 
for 4 of the Galactic globular clusters in the sam- 
ple, where we have normalised the velocity dis- 
persions to their asymptotic values, and the radial 
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0.5 1 1.5 2 2.5 

R/R„ 

Fig. 6. — The curves represent the projected ve- 
locity dispersion profiles normalised to cr and to 
Ro of our models for four representative examples. 

coordinate to the values where the condition x = 1 
(a = do) is met in each of the models, Rq. Mea- 
surements with confidence intervals appear repre- 
sented by different symbols for each GC. We see 
again that <r p (R) profiles consistent with obser- 
vations are easily obtained, while fitting simulta- 
neously surface brightness profiles, to within the 
error bars, where we have also assumed a 20% un- 
certainty in all distance determinations. The four 
clusters included in this figure illustrate the range 
of behaviours found, from NGC 7078 which shows 
a very large variation in a p (R) from the centre to 
the outskirts, implying a large Newtonian inner 
region and only a slight MONDian outer zone, to 
clusters with only small variations in <J p (R) like 
NGC 6171 and NGC 7099. An extreme case is 
NGC 288 with a <r p (R) profile which is also con- 
sistent with observations, but which can not be 
plotted in figure (6), as that low density cluster 
lies in the a < clq regime throughout, and there- 
fore Rq is not defined for it. 

Notice that no external field effect of MOND 
was included in the dynamical modelling per- 
formed, as was also not included in e.g. Sollima 
& Nipoti (2010). That accurate mass and velocity 
dispersion models can be thus constructed, is sug- 
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gestive of a modified gravity formulation where no 
such effect appears. 

Finally, in figure (7) we plot inferred values of 
the total masses and asymptotic values of pro- 
jected velocity dispersions for the 8 globular clus- 
ters we study, with their corresponding confidence 
intervals. The solid line shows the best fit power 
law dependence for this plot, excluding Omega 
Cen. from the fit. Notice that the quantities plot- 
ted here are completely independent of any dy- 
namical modelling, being merely measured veloc- 
ity dispersion values in the outskirts of the clusters 
in question, and total masses as inferred from total 
luminosities and models of their observed stellar 
populations. 

Despite the small ranges in parameters covered 
by the clusters studied, it is clear that the data 
are consistent with the generic modified gravity 
prediction of ao oc Af 1 / 4 , the slope of the solid 
line being statistically determined as 0.31 ± 0.06. 
For comparison, the dashed line gives the best 
slope 1/4 fit, as the numbers above show, a con- 
sistent description of the data. Omega Cen. was 
excluded from the fit, as it is usually acknowl- 
edged as an outlier. The complex stellar popu- 
lations it presents, its abnormally high mass, to- 
gether with the possible presence of an intermedi- 
ary mass black hole in its centre probably indicate 
a complex dynamical formation history. Still, its 
inclusion in the power law fit of figure (7) only 
slightly modifies the resulting slope to 0.32 ± 0.04. 

Figure (7) lends credibility to the interpreta- 
tion of the outer flattening of the observed veloc- 
ity dispersion profiles as evidence for modifications 
in the law of gravity. Under classical gravity, one 
needs to invoke further causal correlations to ac- 
count for the clear trend appearing in figure (7), 
a more contrived scenario than simply the natural 
consequence of a shift towards MONDian dynam- 
ics in globular clusters at scales where a < ao . As a 
comparison, the dotted line in figure (7) shows the 
prediction of eq. (6) for isothermal systems fully in 
the a << ao regime. 

That the GCs treated lye somewhat below the 
prediction of eq.(6) is not surprising, as they are 
by no means structures fully in the a << ao 
regime, most of the mass is partly in the New- 
tonian regime, or in the transition between both, 
hence the asymptotic velocity dispersions slightly 
off the results of eq.(6). 




4 5 6 7 

log(M/Mj 



Fig. 7. — The points with error bars represent 
the values of <7o and total stellar mass for the 8 
globular clusters we study. The solid line gives the 
best fit power law, having a slope of 0.31 ± 0.06, 
excluding Omega Cen. from the fit. If Omega Cen 
is included in the fit, the slope comes to 0.32±0.04. 
The dotted curve shows the prediction of eq.(6) for 
isothermal systems fully in the a << a regime, 
while the dashed line gives the best slope 1/4 fit. 

We acknowledge that a number of the assump- 
tions going into the modelling remain as such, it 
is entirely plausible that some might not strictly 
apply. As already discussed, the case of NGC7078 
could be indicative of a certain degree of orbital ra- 
dial anisotropy appearing towards the outskirts of 
some of these systems. Also, the projected M/L 
values might have some radial variations, or the 
flattening of the projected velocity dispersion pro- 
files might simply reflect the presence of contami- 
nating unbound stars, or the effects of tidal heat- 
ing, as assumed under Newtonian interpretations 
of GC structure e.g. Drukier et al. (2007) or Lane 
et al. (2010). Under Newtonian dynamics, per- 
fectly self-consistent models can be constructed, 
often leading to even better fits than those result- 
ing from modified gravity approaches, e.g. Ibata 
et al. (2011). It is the new results of figure (7), 
which are independent of any gravitational mod- 
elling, which we find most suggestive of the modi- 
fied gravity interpretation, as under the standard 
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Tabic 1: Model parameters for the globular clusters treated. 


Globular Cluster 


log(po) 


o\ 


CO 


r a 


{M/L) M 


{M/L) GC 


NGC 288 


2.15 


0.83 


2.1 


18.0 


2.0 


3.03 ± 1.12 


NGC 1851 


5.65 


6.00 


3.9 


7.0 


2.0 


3.00 ± 1.19 




4. 1 U 


q nn 


o n 
z.u 


o n 
y.u 


1 


o 71 -U 1 nn 


NGC 5139 


3.60 


8.90 


7.4 


15.0 


2.0 


2.68 ±0.98 


NGC 6171 


3.52 


1.20 


2.8 


5.0 


2.5 


3.20 ± 1.20 


NGC 6341 


4.68 


3.50 


3.3 


9.0 


1.63 


2.55 ±0.95 


NGC 7078 


4.98 


5.80 


2.9 


11.0 


1.7 


2.51 ±0.90 


NGC 7099 


4.92 


1.80 


2.3 


8.0 


1.65 


2.60 ±0.90 



po gives the central values of the stellar density used in each model in units of Mqpc~ 3 , while ao, <Ji and r a 
give the parameters of the volumetric velocity dispersion profile used, in units of km/ s and pc, respectively. 
(M/L)m gives the mass-to-luminosity ratios used in each model, and (M/L)gc gives the corresponding 
values reported by McLaughlin & van der Marcl (2005) for each cluster through single stellar population 
modelling using age and metallicity parameters as appropriate for each. 



approach, further mechanisms must be invoked to 
explain the appearance of a " Tully-Fisher" rela- 
tion for GCs. 

We conclude this section with table (1), which 
summarises the observed parameters used and the 
model parameters fitted for each of the 8 galac- 
tic globular clusters in this paper. We note also 
that the models are relatively insensitive to small 
changes in the values of po used, only changes of 
factors <~ 3 and above in this model parameter 
result in significant changes. 

4. Conclusions 

We show that for Galactic globular clusters, 
spherically symmetric equilibrium models can be 
constructed using a modified Newtonian force 
law which reproduces the MOND phenomenol- 
ogy, which naturally satisfy all observational con- 
straints available. 

The resulting models are typically characterised 
by a Newtonian inner region, which smoothly 
transits to a modified gravity outer region on ap- 
proaching the a — a threshold. The resulting in- 



ternal velocity dispersion profiles correspondingly 
transit from an inner radially decaying region to 
an outer flat velocity dispersion one. 

By comparing against careful single stellar pop- 
ulation modelling of the globular clusters studied 
to derive total stellar mass estimates, we show 
that the asymptotic values of the measured ve- 
locity dispersion profiles, a p (R — > oo), and total 
masses for these systems, M, are consistent with 
the generic modified gravity prediction for a scal- 
ing cr^(R oo) oc M. 
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